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The continued-fraction (CF) method was developed 
systematically by Risken and co-workers to solve prob- 
lems of arbitrary fluctuations in nonlinear systems 1 . 
However, this efficient technique is limited to problems 
with a few variables, 2 which in practice means systems 
of noninteracting entities (particles, spins, etc.) In this 
communication, we illustrate how to extend the CF 
method to weakly coupled systems with the problem of 
relaxation in classical spins. 

The Fokker-Planck equation governing the T ^ dy- 
namics of a classical spin s is given by 3 
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where A is the dissipation parameter, B = — (dTt/ds) the 
effective field, and td oc T -1 the free diffusion time. If 
we write this equation as dW/dt = CppW, the dynami- 
cal equation for any average (/) is given in terms of the 
adjoint Fokker-Planck operator by d (/) /dt = (4p/)- 

However, the underlying Landau-Lifshitz damping 
term —As A (s A B) is nonlinear in s and couples the 
equations for the moments. To handle the infinite sys- 
tems of coupled equations it is convenient to introduce 
the spherical harmonics X™ = e lmv P™ (s z ) , and take 
advantage of their recurrence relations. Writing £ FP in 
terms of the rotation operator J — —is A (d/ds) 


4p = -jB- [J + X{sAJ)] - J 2 


(1) 


we get the equations for the X™ in the form of a recur- 
rence relation that can be solved by CF methods. 

When the spins interact, the structure of £ FP is the 
same with the field created by the other spins bd included 
in B. This allows to separate the free evolution and inter- 


action terms, £ 
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bd ■ Jll, where J hL is the 

operator multiplying B in Eq. (I). For weak coupling, 
we can expand the averages in the coupling parameter 
(/) = (/)o + + " " •> getting a hierarchy of equations 
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needs (X™Y£)., i = 0,1, which could be obtained writ- 
ing the above equations also for / = X" 1 Y£. This leads 

to another complication, since bd ■ Jll f then involves 
products of three spherical harmonics. However, Y^Y^ 
can be expanded as a linear combination of the Y^'s, 
restoring a structure solvable by CF methods. 

These manipulations result in a complicated algorithm, 
so it is convenient to check it out against some solvable 
case. Fortunately, for freely rotating spins the problem 
reduces to that of electric dipoles, where the perturbative 
calculation can be done analytically, yielding 4 
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where ao(w) and a (oc T _1 ) are the dynamic and static 
polarisabilities of noninteracting dipoles, g is the dipole 
density, and 1Z and S are certain "lattice sums" of the 
interaction tensor. Figure 1 shows that our perturbative 
CF approach gives accurately the analytical results. 



FIG. 1. Cole-Cole plots of the dielectric polarisability 
(imaginary vs. real part) for various couplings q. Lines: 
Zwanzig formula; crosses: perturbative CF results. 

Therefore, this approach extends the applicability of 
the efficient CF techniques to (weakly) interacting sys- 
tems, and allows one to handle problems involving driven 
systems, evolving in nonlinear potentials, coupled, and 
subjected to dissipation and fluctuations 
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To illustrate our procedure, let us denote by X™ and 
Yp 9 the spherical harmonics of two of the spins, and con- 
sider the equation for X™ (i.e., / = X™ above). Since 
the field created by the second spin depends on Y^, one 


* wzar . unizar . es/ acad/f ac/cie/ cond-mat/T/ theory . html 

1 H. Risken, The Fokker-Planck Equation (Springer, 1989). 

2 P. Hanggi et al, Rev. Mod. Phys., 62, 251 (1990) 

3 W. F. Brown, Jr., Phys. Rev. 130, 1677 (1963). 

4 R. Zwanzig, J. Chem. Phys. 38, 2766 (1963) 


